Protocols for characterising quantum transport through nano-structures 
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In this work, we have analysed the exact closed-form solutions for transport quantities through 
a mesoscopic region which may be characterised by a polynomial functional of resonant transmis- 
sion functions. These are then utilized to develop considerably improved protocols for parameters 
relevant for quantum transport through molecular junctions and quantum dots. The protocols are 
shown to be experimentally feasible and should yield the parameters at much higher resolution than 
the previously proposed ones. 
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The single particle scattering approach, pioneered by 
Landauer [J, Q , and extended by Biittiker Q has been 
the most extensively employed framework for investi- 
gating quantum transport through nano-structures. A 
steady-state non-equilibrium problem is mapped to a 
scattering problem in this approach. Realising the impor- 
tance and applicability of such an approach, a Landauer 
like formula was derived by Meir and Wingreen [|| for an 
interacting mesoscopic region coupled to non-interacting 
leads with coupling strengths of IY and Tr to the left 
and right lead respectively. For IY = AFr, the current 
through the region flowing into one of the leads may be 
expressed as an integral transform, namely, 



a non-interacting region, eq (JTJ) reduces to the Landauer 
formula in which case the L(e) is simply the transmission 
function of the device under consideration. 

For a specific form of L(e), such as a resonant trans- 
mission function (RTF) (a Lorentzian) , an exact solution 
of cq exists Thus, if the L(e) can be represented 

as a polynomial functional of RTFs, such as a linear su- 
perposition of multiple (Lorentzian) resonances (at the 
lowest order) 
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in terms of the local properties of the region (L(e)) by 
a kernel K(e; V) over the real line. The kernel is given 
by K(e;V) = f(e - eV/2;T L ) - f(e + eV/2;T R ), with 
f(x; T) = ( e x / k BT _|_ 1 ) l ag t i ie Fermi-Dirac distribu- 
tion function; Tl and Tr are the temperatures of the 
left and right reservoirs respectively, and V, distributed 
symmetrically across the two electrodes as the voltage 
bias. The quantity representing the device region is 
L(e) and for interacting systems, is represented in terms 
of full non-equilibrium Greens functions (retarded (G r ) 
or/and advanced (G a )) of the interacting device region, 
in close resemblance with the Landauer formula as 
L(e) = ^Im[tr{rG r }], where, T = r L r R /(r L + T R ). 
Even if Tl 7^ ATr, it was shown by Ness et. al. [f| that the 
current can be written as eq ([T]) with renormalized cou- 
pling to the electrodes. Such a Landauer-like approach is 
suitable for any mean-field based method, for example, 
density-functional-based techniques (DFT, TDDFT) or 
interactions treated at the Hartrcc-Fock level For 



the exact solution of eq ([T]) can be obtained (with the 
use of partial fractions for higher orders). Such a form 
has indeed been found for a variety of nano-systems such 
as molecular junctions, quantum dots and quantum point 
contacts [a-tl2| , where the discrete level e r gets broadened 
due to the coupling (7,.) to the macroscopic leads, and 
A r is an unknown parameter that depends on a number 
of factors like the number of conducting molecules in the 
junction or the coupling with the electrodes [l3|. More 
importantly, it is not restricted to systems in the ballis- 
tic regime and is applicable even for interacting meso- 
scopic systems under a broad range of experimen tally 
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relevant conditions [H [^, UH, [l4| . As mentioned in 
such a Breit-Wigner (BW) type resonant conductance 
formula is relevant for an interacting system with a non- 
degenerate ground state like in semiconductor nanostruc- 
tures or ultra small metallic systems. The positions and 
intrinsic widths of the BW type resonances are deter- 
mined by the many body states of the interacting sys- 
tem. The magnitudes of the temperature (ksT), bias eV 
and coupling to the electrodes should be much smaller 
than the resonant energy so that only a single transi- 
tion from the N electron ground state (GS) to N + 1 
electron GS is allowed. At a finite voltage these reso- 
nances may get shifted [l3l4l~5j relative to the zero bias 
position. This shift manifests itself differently for a sym- 
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metric or asymmetric junction. With minor change of 
e r — > e r + rjV, asymmetric couplings to the electrodes 
may also be incorporated with 77 = being the symmet- 
ric case. Naturally, strongly interacting systems such as 
those where Kondo physics is important exhibit slow log- 
arithmically decaying tails in L(e), and hence cannot be 
captured within such an approach [l6| , since Lorentzians 
have a algebraically decaying tail structure. 

In this work, we have analysed the exact solutions of 
eq ([T]) with an L(e) given by cq ([2]) and discuss their gen- 
eral validity. Further, we utilize them in designing sub- 
stantially improved protocols for finding the parameters 
of the L(e), specifically for quantum dots and molecular 
junctions. These protocols are shown to be experimen- 
tally feasible and if implemented, will yield the parame- 
ters with much higher resolution than the existing ones. 
An understanding of the non-linear regime, both in terms 
of voltage and thermal bias is important. Such an insight 
is most easily developed through closed form analytical 
expressions. The asymptotic properties of exact solutions 
are the most straightforward route to obtaining such ex- 
pressions, thus emphasising the utility of exact solutions. 

Substituting eq ((2|) with A r = 7^ in eq dJ), and trans- 
forming it to a contour integration Q , we get the follow- 
ing expression for the current I(V,Tl,Tr). 
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where ^f(z) is the digamma function jl7| and z^m^ = 
1/2 + [j r - i(e r =F eV/2)]/(2nk B T L/R ). Using the current 
expression thus obtained, we can simply write down the 
differential conductance G(V,Tl,T r ) = dl/dV and dif- 
ferential thermopower S(V, T L , T R ) = dAV/dAT. These 
are given by 
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where ^'(z) is the trigamma function [l7j and 
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An expression for thermopower in terms of trigamma 
functions has been obtained earlier in the linear response 
regime \l§L [l9j . The expression derived here is exact 
and hence represents a generalization of that result to 
all regimes. 

It is worth considering the general structure of the 
equations above. The thermal energy, k B T L / R sets 
the reference scale, since the final expression is a func- 
tion Only Of £ L /R,r = ^L/R,r/ {^k B T L / R ), 7L/_R,r = 

j r /{2nk B T L/R ) and V = eV/(2nk B T L/R ). The width 
of the resonance appears in the real part of the argu- 
ment, while the peak energy and the bias voltage appear 
in the imaginary part of z^/ R . Since many of the proper- 



ties of the digamma function [20( depend on the real and 
imaginary parts separately, we can classify the param- 
eter space into the following regions: narrow resonance 
(7V <C 1), finite width resonance (tv ~ 1) and broad 
resonance (tv ^ 1). 

a. Narrow resonance: The narrow resonance regime 
is characterized by having the width of the resonance 
being much smaller than the thermal energy scale (tv <C 
I). In the expression given by equation if j r <C 1, 
then using the identity, Im\l/(l/2 + iy) = (7rtanh7r?/)/2 
[13, we get 
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where a L / R r = (e r ^feV/2)/k B T L / R . A similar result has 
been obtained before, through Keldysh formalism [2l[ for 
resonant transmission through one dimensional systems. 
The differential conductance in this regime can be ob- 
tained by calculating G = dl/dV and is given by, 
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which yields conductance oscillations function of 

source-drain bias (when eV = ±2e r ) or as a function 
of gate voltage (which tunes e r ) at zero bias [ll|, [HI- It 
is well known that the positions of resonances in the zero 
bias conductance versus gate voltage curve yields values 
of the resonance energies. The interpretations of these 
oscillations as being due to resonant transmission or due 
to Coulomb blockade rests on the dependence of the en- 
ergy level spacing AE r = e r +\ — e r on the bias. If the 
spacing increases with increasing bias, then the energies 
are single-particle energies, while for constant spacing, 
the levels are many-particle levels that include charging 
energy (T2T [T3|. 

If the individual resonance peaks are separated by 
energies far greater than the thermal energy scales 
{k B T L / R ) and the widths (7,.), then the resonance clos- 
est to the chemical potential would be the major con- 
tributor to the current, and hence a single resonance TF 
may be assumed with width 7 and peaked at eo- Such 
a situation may be realized in quantum dots by reach- 
ing sufficiently low temperature. In such a case, we can 
obtain the well known expression for thermo voltage, Vth, 
as, Vth = -Sideai®, where, S lde ai = eo/(eT) is the ther- 
mopower of an ideal quantum dot characterised by a zero 
width (5- function) TF, T is the average temperature 
given by T = (T R + T l )/2 and 6 = T L - T R is the ther- 
mal bias. 

In a recent work by Mani et. al. [23[, a protocol for 
obtaining the width, 7, of an RTF was proposed through 
the measurement of a quantity termed thermopower off- 
set (F) defined as F = (S ide ai - S)/S ldea i- This quantity 
was shown, employing numerical calculations, to be a 
simple polynomial function oi^/k B T (when eo — > 0), i.e. 
F = A(j/k B T) + B(-//k B T) 2 with A and B being con- 
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stants. By an experimental measurement of F, the above 
equation can be inverted to find ^y/kgT. 

Since this protocol relies on an accurate measurement 
of S in the limit eo — > 0, where S would itself be vanish- 
ingly small, F would be a difficult quantity to measure 
with high resolution. Hence we propose an alternative 
protocol for finding 7, which does not require tuning of 
the resonant level to zero. We first state that the ther- 
mopower of a quantum dot in the linear response regime 
and in the limit eo — > may be represented as, 
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where A,B and C are pure constants [2J], given by 
A = 3C(2), C = 14C(3) and B = 2C (((z) is the Reimann 
zeta function), which implies that the linear term in the 
offset expression has the coefficient 7C(3)/37r£(2) ~ 0.54, 
which matches with the value obtained by Mani et al [23| 
through numerical fitting, and also shows that these co- 
efficients are indeed pure numbers. For eo 7^ 0, we obtain 
a general expression for the thermopower [24|, given by, 
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FIG. 1. Main panel: The relative discrepancy in ther- 
mopower, AS/S, denned as (Sexact — S pr otocoi)/ S exac t (com- 
puted using equations [S] and [5} is plotted as a function of 
7/fcsT for non-zero eo = 47rfcsT (blue solid line) and 2-KksT 
(red dashed line) at T = 3 K and 6 = 0.1 K. Inset: The same 
quantity is calculated in the nonlinear regime at T = 3 K 
and 9 = 3 K [2(|. The relative discrepancy reduces to 1-2% 
in the nonlinear regime thus offering an improved standard. 
The schematic represents the relevance of this protocol for 
quantum dot devices. 
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The quantity is the general ideal thermopower (in 
units of fcs/e), that reduces to eo/eT in the linear re- 
sponse regime. The quantities Co and C\ are functions 
purely of eo, Tl and Tr and can be expressed in terms 
of polygamma functions. These may be easily evaluated 
either using the series expansions given in the SM [24[ 
or through technical software like MATLAB [25| . Thus 
the protocol simply consists of measuring thermopower 
at a given eo, which can be chosen such that high reso- 
lution is achieved; finding the coefficients Fx, F , Cq and 
C\ using the expressions given in SM [24j], and inverting 
equation © or equation to get the resonance width 
7- 

We have bench marked the expression used for this 
protocol in Fig. ([IJ, where we show the difference be- 
tween the protocol expression equation © and the ex- 
act result equation (|5|) as a function of the scaled width 
7/fcgT for various resonant level positions. It is seen 
that the discrepancy, AS/S, between the exact (equa- 
tion and the protocol expression (equation ©) is 
less than 4% for an eo as large as 27rfcsT. In fact an 
analysis of thermopower tells us that the discrepancy, 
AS reduces dramatically in the nonlinear regime thus 
offering an even better protocol. This is seen in the inset 
of fig. ((l) where AS is seen to decrease in the presence 
of a finite thermal bias (as compared to the main panel) . 
We have seen that the discrepancy between the protocol 
and exact thermopower offset near eo — > can be as large 
as 20%, and is hence less reliable (see Figure (4) in SM 
M). 



b. Broad Resonance: The broad resonance trans- 
mission defined by 7 ^> UbTl, ksTn is most appropriate 
for molecular junctions where the HOMO and LUMO 
levels are broadened due to the coupling with the leads, 
and the width of these levels could easily be of the or- 
der of eV, which is far greater than the thermal energy 
scales. Further, since the HOMO-LUMO level separation 
is much greater compared to the width or the thermal en- 
ergy scale, a single resonance £(e) can again be assumed. 
It is easy to see that in equation if Tl/r S> 1, then 
l^i/i?! ^> 1 irrespective of the values of the eo or the bias, 
since the latter are present in the imaginary part, while 
7 is in the real part of Zl/r. This allows us to use the 
asymptotic form of digamma function (l7| for large z, 
which is ^(z) ^ \nz — l/z;z— > 00, |arg(z)| < ir. Hence, 
we get for / in units of 2e/h, 



I = 



tan 



lr~ K l 



T L X r 



tan 



A r -| 

7r 



( 7r 2 + A 2 J 2 (7 2 + A? + ) 2 



(10) 



where, A r ± = e r ± eV/2. As the above arguments show, 
this equation is valid for arbitrary values of the reso- 
nance position e or the voltage bias V as long as 7 1 
is satisfied implying that it is sufficient for the thermal 
energy scales to be much smaller compared to the reso- 
nance width. Although the above arguments seems to be 
based upon a broad resonance condition, a subtle point 
to note is that this form may be applied for arbitrary 
widths (pf/hsT) if the magnitude of the level position 
(measured from the chemical potential) is large compared 
to max(eV/2, ksT), thus making it useful for molecular 
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junctions. The result obtained here represents a general- 
ization of an expression obtained by Stafford [H|,[27]], (for 
RTFs in the wide band approximation) at T = to fi- 
nite temperatures and a finite thermal bias. We can now 
clearly see through eq (fTO)) , the emergence of a temper- 
ature controlled current rectifier. This rectification cur- 
rent, defined as AI = (I(V,T L ,T R ) + I(-V,T L ,T R ))/2 
is given by, 
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This rectification current is experimentally measurable 
(~ nA) (Figure (3) in SM (28|) even at a temperature and 
thermal bias of 50 K each. This motivates us to design 
a protocol for predicting the resonant transmission func- 
tion parameters through AI at zero voltage bias (which 
is in fact the thermocurrent, AI t h) in conjunction with 
the zero bias conductance. This protocol involving the 
thermocurrent (AI t h) shall be discussed below. 



The existing protocol for finding the resonant level in 
molecular junctions is implemented through transition 
voltage spectroscopy (TVS) p9M33| . The basis for this 
protocol is the existence of a minimum, in the In (I /V 2 ) 
versus l/V curve. This minimum is useful because, it oc- 
curs at a voltage that is much smaller and accessible than 
the resonance condition voltage (V = 2eo). The interpre- 
tation of this voltage minimum as V™ nv = 2e / -\/3 relies 
on a result [27| obtained in the limit when 7/eo "C 1 , and 
hence the information on the width of the resonance is 
completely lost. We have obtained a result for the TVS 
minimum that is valid for a regime where 7 /en may be 
significant (e. g. amine linked junctions) [34L [35j, which 
is, 
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A comparison between the V m obtained and the V^ nv 
with the exact result [28j as obtained by numerically find- 
ing the minimum is shown in fig. ([2]). It is seen that 
the deviation from the conventional TVS minimum oc- 
curs even at very small widths (7/eo ~ 0.05), while the 
agreement with the exact result is excellent. Thus the 
experimentally measured V m would contain information 
about both e and 7 and hence by itself, cannot be used 
to find the level position and the width. A second equa- 
tion is needed that relates an experimentally measurable 
quantity to eo and 7. We state that such a quantity is 
the ratio Rtvs of AIth to the zero bias conductance G. 
The latter (obtained either by using equation (j4|) or using 
equation (fTU|) . for a single resonance with Tl = T R = T 
is given by 
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FIG. 2. Comparison of TVS minimum, Vm, obtained us- 
ing the protocol equation (|12[) (red dashed line), exact so- 
lution [23] (black solid line) and the conventional [23 Vm™ = 
2eo/%/3 (blue dashed line). The schematic in the inset repre- 
sents a molecular junction. 



and hence the ratio is given by: 
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Equations (|T2"j) and (fl"4"]l can be easily combined to get 
both 7 and eo- So if we define X — v / 3eV m /2 and 
Y = 3ei?TVs/(47r 2 T8), which are experimentally mea- 
surable, then 7 and e may be obtained through simple 
expressions involving X and Y (28j . 

c. Finite width: In the final part, we provide ex- 
pressions for the I — V characteristic in the finite width 
regime (7 ~ ksT), which has hitherto been analyti- 
cally inaccessible. We can utilize the recurrence relations 
and the multiplication formula for the digamma function 
[13 to get the following result. If T L = T R = T and 
7,. = mtksT] n = 1,2,..., then the current / in units of 
2e/h is given by, for n = 2m + 1 
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The case for 



2m can also be 



where, v L/R = - — ; 
similarly derived [361 ] . 

Finding the L(e) for systems where electron-electron 
and electron-phonon interactions play a dominant role, 
is of course a major bottleneck and is currently a ma- 
jor topic of research [5L I37I - I42I ] . Nevertheless, it is clear 
from our work that expressing the local properties of the 
interacting region (L(e)) in terms of RTFs allows the 
utilization of exact results. Subsequently, we have pro- 
posed substantially improved protocols that can be im- 
plemented experimentally for finding transmission func- 
tion parameters in quantum dots and molecular junc- 
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tions. Our results being based on exact solutions also 
provide analytical insight into the difficult nonlinear 
regime, and provide a unified platform for the analy- 
sis of simulations and experiments in quantum transport 
through nanostructures. 
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